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Upon employing a natural time window of fixed length sliding through a time series, an explicit 
interrelation between the variability /3 of the variance «i(= (x 2 ) — (x) 2 ) °f natural time x an d events' 
correlations is obtained. In addition, we investigate the application of the fluctuation theorem, 
which is a general result for systems far from equilibrium, to the variability /?. We consider for 
example, major earthquakes that are nonequilibrium critical phenomena. We find that four (out of 
five) mainshocks in California during 1979-2003 were preceded by j3 minima lower than the relative 
thresholds deduced from the fluctuation theorem, thus signalling an impending major event. 
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Entropy production is a measure of the irreversibility 
of a thermodynamic process: the difficulty, even impos- 
sibility, of reversing the observed often macroscopic be- 
havior of a system that exchanges heat or matter with 
a complex environment (e.g., see Ref.JH and references 
therein). The breakage of time reversal symmetry as- 
sociated with thermodynamic irreversibility has focused 
enormous discussion for more than a century. Despite 
of such concerns, however, the concept of entropy gener- 
ation in the thermodynamics of large systems has been 
applied widely. From microscopic point of view, efforts 
towards understanding the nature of the entropy and its 
production -mainly focused on the one way character of 
the second law- have been attempted. They modelled the 
microscopic evolution of a system and its environment in 
the frame of stochastic dynamics and stochastic ther- 
modynamics but interpretations based on deter- 
ministic dynamics (e.g., see Ref.@) were also forwarded. 

An intense interest towards the latter interpretations 
has been renewed when Evans, Cohen and Morris in 1993 
considered the fluctuations of the entropy production 
rate in a shearing fluid, and proposed the so-called fluc- 
tuation relation or the first fluctuation theorem Q . This 
is considered [1| to represent a general result concerning 
systems arbitrarily far from equilibrium. The proof of the 
fluctuation Q and related theorems shows how irre- 
versible macroscopic behavior arises from time reversible 
microscopic equations of motion. The two theoretical 
results that illustrate this clearly are the second law in- 
equality fill ] and the very recent mechanical proof (l2l | 
of Clausius' inequality without the prior assumption of 
the second "law" of thermodynamics. These two results 
have been obtained without treating the nonequilibrium 
entropy, but used instead a quantity termed dissipation 
function first defined [l3| in 2000. On the basis of this 
function, being a path function and not a state function, 
the relaxation of a system to equilibrium, which is inher- 
ently a nonequilibrium process, can be quantified [3]. 

It has been emphasized in Ref. @ that, unlike linear 
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irreversible thermodynamics, the fluctuation and related 
theorems are exact for systems of arbitrary size as well 
as for systems arbitrarily near to, or far from equilib- 
rium, as mentioned. This is why we shall employ here 
the fluctuation theorem for the purpose of the present 
study. 

This theorem 0, [H|-[l9| gives a general formula for 
the probability ratio that in a thermostatcd dissipativc 
system, the time average entropy production takes a 
value A to minus the value —A, 

Pr{t t /k B = A)/Pr{t t /k B = -A) = cxp[At] (1) 

from which it is obvious that as the averaging time or sys- 
tem size increases, it becomes exponentially likely that 
the entropy production will be positive. The theorem 
was initially proposed Q for nonequilibrium steady states 
that are thermostatcd in such a way that the total energy 
of the system is constant. Subsequently, it was shown 
fl8l [l9j that this theorem can be proved for sufficiently 
chaotic, iso-energetic nonequilibrium systems using the 
Sinai-Ruelle-Bowen measure, as well as for purely Hamil- 
tonian systems with or without applied dissipative fields 
[20l | and for a wide class of stochastic nonequilirium sys- 
tems nnm. 

It is one of the two basic aims of this paper to investi- 
gate for the first time the application of the fluctuation 
theorem to the case of earthquakes which may be consid- 
ered (e.g. [23l |2~H) as nonequilibrium critical phenomena 
(the mainshock being the new phase) . They exhibit com- 
plex correlations in time, space and magnitude M which 
have been recently studied by several workers (e.g., see 
Refs. [25rl29l ] ) . In particular, the present investigation 
will be made by applying the fluctuation theorem to the 
order parameter fluctuations that result from the anal- 
ysis of the time series in a new time domain termed [30l| 
natural time x- This is so, because natural time analy- 
sis allows us to identify [3l| when a complex system ap- 
proaches a critical point (for a review see Ref. (HI) and 
in addition enables the introduction of an order parame- 
ter for seismicity. The present study has been motivated 
by the following two findings related to the variability B 
(defined below) of the order parameter of seismicity [331 ] : 
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First, it captures the events' correlations, as shown here 
(see Appendix) , which constitutes the other basic aim of 
this paper. Second, the quantity /3 exhibits characteristic 
minima }34j before the occurrence of major events. 

In a time series comprising TV earthquakes, the nat- 
ural time Xk = k/N serves as an index for the occur- 
rence of the fc-th earthquake. In natural time analy- 
sis the pair (xkiQk) is studied, where Qk is the energy 
released during the fc-th earthquake of magnitude Mk- 
One may alternatively study the pair (XkiPk)i where 
Pk = Qfc/En=i Qn is the normalized energy released 
during the fc-th earthquake, and Qk -and hence pk- is 
estimated through the relation [H] Q k oc 10 1,SM *. The 
variance Ki(= (y 2 ) — (x) 2 ) of x weighted for pk, is given 
by P, H, m 



N 



N 



«x = ^2pk(xkf - ^2?kXk 



(2) 



k=l 



v.fc=l 



This quantity, as shown in Ref. [33[, can be also consid- 
ered as an order parameter for seismicity. 

The fluctuations of Ki, are studied by applying the 
following procedure [32[ • Taking an excerpt of a seismic 
catalog comprising W{> 100) successive events, we start 



from the first EQ and calculate the first 35 K\ values for 
6 to 40 consecutive EQs. Then we proceed to the second 
EQ, and calculate again 35 values of K\ from the 7-th to 
the 41-st event. Thus, scanning event by event the whole 
excerpt of W earthquakes, we calculate the average value 
//(ki) and the standard deviation cr(/ti) of the K\ values. 
The quantity 



= o-(aci)//z(ki) 



(3) 



is defined [33 as the variability /3 of Ki for this excerpt 
of length W. In some occasions, as in the present case, 
it is of prominent importance to know what happens to 
the /3 value until just before the occurrence of each EQ, 
ei, in the seismic catalog. We then calculate first the 
Ki values using the previous 1=6 to 40 consecutive EQs. 
These 35 K\ values are associated with the EQ e,, but 
we clarify that EQ ej has not been employed for their 
calculation. The j3 value -corresponding to the EQ ej- for 
a natural time window length W is computed using all 
the (35 x W) K\ values associated with the EQs ei-w+i 
to ej. The resulting value is denoted by (3yy, where the 
subscript W shows the natural time window length, and 
the corresponding minimum is designated by (3w,min ■ 
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It is shown that the quantity j3 when using I consecutive events is interrelated with the event's correlations through 



-E 



all pairs 



(x) 



M I 



(x> 



M , 



Co\{p j: p n 



E 



all pairs 



(j ,r )2 cov(p j ,p m ) 



K hM + Eall pairs ^^Covipj , Pm ) 



(4) 



where (x)m an d Ki t M correspond to the average value of x an( i K ii respectively, obtained when substituting for pk 
the average -within an excerpt of W events- values [ik of pk ; the symbol Cov(pj,p m ) stands for covariance, i.e., the 
average value of (j>j — p,j)(p m — Mm) within the excerpt of W events. The details of the derivation of Eq.(j4]) are given 
in the Appendix. 



The selection of the W value used for the purpose of 
our study is of crucial importance. It is taken equal 
to the number of the events that would occur in a few 
months, or so, in view of the following: Low frequency 
(< 1 Hz) electric signals, termed Seismic Electric Sig- 
nals (SES), appear before earthquakes l39l l40j . They are 
emitted from the future focal region [4JJ (see also Ref. 
[42| ) when in the focal region the stress reaches a crit- 
ical value <7 cr , and then a cooperative orientation of the 
electric dipolcs occurs. This leads to the emission of a 
transient electric signal that constitutes an SES. Several 
suc h sig nals within a short time are termed SES activity 
[HI, l37l. |42T |43| . For example, the three lower channels 
in FigOTb) show three SES activities that preceded ma- 
jor earthquakes in western, southwestern and southern 
Greece, respectively, as depicted in the map of FigfJTa). 
Furthermore, for the sake of comparison, the upper chan- 
nel in Figfljb) shows a very recent SES activity initiated 
on 8 January 2013 at a station labelled LAM in Fig. Ufa) 



in central Greece. The following important fact has just 
been identified 441: At the initiation of an SES activ- 



ity, which usually occurs a few months (with an upper 
limit of around 5 months) before a major EQ, a clearly 
detectable change in seismicity appears, manifested by a 
minimum Pw,min in the fluctuations of the order param- 
eter of seismicity. Hence, in the case that geoelectrical 
data are lacking, once we identify the date of pw,min (by 
analyzing solely seismic data) this reveals also the date 
of an SES activity that would have been recorded. 

Along these lines, Table U shows the dates of 
the minima /3w,min °f seismicity before major main- 
shocks in California during the 25 year period 1 
January 1979 to 1 January 2004. We used the 
United States Geological Survey Northern California 
Seismic Network catalog available from the Northern 
California Earthquake Data Center, at the http ad- 
dress: www .ncedc . org/ncedc/ catalog-search. hmtl, 
hereafter called NCEDC. The seismic moment Mo, which 
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FIG. 1: (color online) (a) Major earthquakes in Greece on 8 June 2008 (green, magnitude M w =6.4), 14 February 2008 (red, 
M w =6.9 and 6.4) and 8 January 2006 (blue, M w =6.7). (b) Their preceding SES activities recorded at Pirgos (PIR) measuring 
station located in western Greece are shown (with the corresponding color) in the lower three channels. Furthermore, an SES 
activity initiated recently on 8 January 2013 at a station in central Greece labelled LAM in (a) is depicted in the upper channel 
of (b). 



TABLE I: The minimum DFA exponent a m in along with the values of the minima observed for the variability j3 together with 
the dates of their observation in parentheses before all major EQs in California with M > 7.0 within N 3 l$ W\^[l during the 
period 1979-2003. The M6.9 Northridge earthquake is also added in italics. The lead time At for each case, estimated from 
the difference in the dates between the EQ occurrence and the appearance of p3oo,min is shown in the last column. The values 



min and a min are taken from Ref. [34| 


EQ Date 


EQ Name 


M 


0300, min 


P200, min 


ttmin 


At (months) 




epicenter 




(date) 


(date) 


(date) 




1980-11-08 


Eureka 


7.2 


0.444 


0.432 


0.445 


« 3 




N41.08°W124.62° 




(1980-08-01) 


(1980-06-28) 


(1980-08-01) 




1989-10-18 


Loma Prieta 


7.0 












N37.04°W121.88° 












1992-06-28 


Landers 


7.4 


0.378 


0.377 


0.383 


< 5 




N34.19°W116.46° 




(1992-01-28) 


(1992-01-03) 


(1992-02-02) 




1994-01-17 


Northridge 


6.9 


0.459 


0.324 


0.431 


w 2 




N34.23°W1 18.55° 




(1993-11-14) 


(1993-10-18) 


(1993-11-14) 




1994-09-01 


Mendocino 


7.0 


0.472 


0.474 


0.458 


« 1 




N40.41°W126.30° 




(1994-08-01) 


(1994-07-11) 


(1994-08-09) 




1999-10-16 


Hector Mine 


7.0 


0.444 


0.432 


0.422 


« 5 




N34.60°W116.34° 




(1999-05-14) 


(1999-05-14) 


(1999-05-15) 




Fluctuation theorem 


and natural time analysis 




0.46 


0.45 







is proportional to the energy release during an earth- 
quake and hence to the quantity Qk used in natural time 
analysis, is calculated |32| from the relation log 10 (Mo) = 
1.5Af+const, where the earthquake magnitudes reported 
in this catalog are labelled with M. The earthquakes 
with M > 2.5 reported by NCEDC, within the area 
Nff^WiH;! have been considered. We have on average 
~ 10 2 EQs per month since 31832 earthquakes occurred 



for the 25 year period from 1 January 1979 to 1 January 
2004. Thus, we adopted natural time window lengths 
W =200 and W =300. 

The results of this analysis are depicted in Fig. [^a),(b) 
where we plot the variability /3 (in red for W = 200 and 
in blue for W = 300) versus the conventional time for 
the periods (a) 1 January 1979 to 1 January 1990 and 
(b) 1 January 1990 to 1 January 2004. An inspection of 
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FIG. 2: (color online) The variability f3 of k\ (left scale) plotted versus the conventional time for a natural time window of 
length W =200 events (red) and W =300 events (blue) during the period: (a) 1 January 1979 to 1 January 1990 and (b) 1 
January 1990 to 1 January 2004. The earthquakes with M > 6.5 (right scale) are shown with vertical bars ending at solid 
circles. 
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FIG. 3: Application of Eq.Q for the experimentally deter- 
mined Pw for W = 200 (red) and 300 (blue). The threshold 
values /3200.0 = 0.45 and /33oo,o = 0.46 are deduced from the 
maximization of the linear correlation coefficient (Pearson's) 
r, thus pointing]!^] to optimal linearity. 



these results lead to the Pw.min values inserted in Ta- 
ble ID In five out of the six mainshocks we find values of 
/?300,mm and (3200,min that appear 1 to 5 months before 
mainshocks. In these five cases ,8200, min varies between 
0.324 to 0.474 and ft 00 , mm between 6.378 and 0.472. We 
note that the key criterion to distinguish the true pre- 
cursory Pw.min from the non precursory ones is the fol- 
lowing |34j : The minimum should be followed (before the 
occurrence of the mainshock) by a period during which 



the exponent a of the Detrended Fluctuation Analysis 
(DFA) [HI -calculated for a length W=300 events in the 
magnitude time series- reaches a minimum a mm slightly 
smaller than 0.5 (thus, indicating anticorrelated behav- 
ior, but close to random) and then ^200 > ^300- This in- 
equality means that when the system approaches closer 
to the critical point -which is the case when considering 
1^=200 events compared to 14 7 =300 events- the fluctua- 
tions of the order parameter become more intense. 

We now proceed to the investigation of Eq.([T]) in the 
case of scismicity and analyze the statistical distribution 
of the experimentally determined @w for W = 200 or 300, 
which is clearly path depended. The quantity (3w can be 
considered as an entropic measure (see Appendix), but 
its sign is by definition always positive. Thus, in order 
to apply Eq. (HJ , we need to define a threshold value /3w,o 
above which the entropy production may be considered 
positive whereas when below negative. For this reason, 
we employ the relation 



Pr {Pw ~ M 
Pr (/3w,o - Pw) 



exp [r' 08 



w 



MY. 



0)] 



(5) 



which results from Eq.([T]) when considering A = /3w — 
fiw.o an d experimentally determine Pr {Pw ~ Pw.o) by 
using bins of width A(3w =0.01. Figure [3] depicts the 
natural logarithm of the left hand side of Eq.© as a 
function of (0 W - p w ,o) for W = 200 and W = 300. 
In each case, the threshold @w,o is the one that maxi- 
mizes the linear correlation coefficient (Pearson's) r, thus 
pointing[45| to optimal linearity. We find the threshold 
values of ,8200,0 = 0.45 and ,8300,0 = 0.46. Moreover, the 
relative 'time-scale' t', which corresponds to the slope of 
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Fig El lies in the range 40.6 to 48.3, which is comparable 
with a scale of the order of I =40 sequential events used 
in the calculation of j3w- 

Let us now compare the aforementioned threshold val- 
ues #200,0 = 0.45 and #300,0 = 0.46 with the (3w,mm val- 
ues identified before each mainshock in Table Q] We find 
that except one mainshock, i.e., the Mendocino EQ in 
1994, the other four mainshocks (including the strongest 
in Table H]) led to #200, mm and /3soo,min values that are 
lower than ,8200,0 and #300,0 j respectively. 

Thus, in summary, it may be considered that once the 
natural time analysis leads to an identification of precur- 
sory minima (i.e., # 2 oo,mm and /?300,min) that are lower 
compared to the threshold #iy,o values determined from 
the combined use of natural time analysis with the fluc- 
tuation theorem, a forthcoming major EQ is likely to 
occur. 



Appendix: Interrelation of the variability f3 with 
correlations when a (natural) time window of fixed 
length is sliding through a time series 

Here, we focus on the mean value [i = of 
Ki and the corresponding standard deviation a = 



1. The mean value /i = £(ki) of aci 

In a window of length / starting at k = fco, the quanti- 
ties 



Qk„+j-i 



X/m=l Qko+m— 1 



1,2,. 



(A.l) 



representing the normalized energy are obtained, which 
satisfy the necessary conditions 



Pj(k ) > 0, 



i=i 



(A.2) 



(A.3) 



to be considered as p oint probabilities. We can then 
define as usual [10, Elf the moments of the natural time 
Xj = j/l as (x q ) = J2 l j=i(j/ l ) q Pj(k ) and hence 



£ [ki — £(ki)] when a (natural time) window of length 
I is sliding through a time series of Qk > 0, k = 
1,2,... W. Once these quantities have been evaluated, 
the variability # of k\ can then be estimated by # = a / 1 \i. 



J 



k iOo) = 



/ / . \ 2 
j 



Pj(k ) 



(A.4) 



Note that n\ is a non-linear functional of {pj}- 



Let us consider the average value fij of pj obtained when the (natural time) window of length I slides through a 
time series of Qk > 0, k = 1, 2, . . . W, i.e., we have[48j 



W-l+i 



w-1+1 



W - 



Qk +j-l 



k = l " ' fe =l Sm=lfto+m-l 

It is obvious that the definition of Eq. (IA.5l) is consistent with Eq. (|A.3|) . thus we have 

I> = 1. 

i=i 



(A.5) 



(A.6) 



Similarly for the second order moments of pj, one can estimate [481] the variance of pj by 



Vax(pj) = £ (pj-fJ-jY 



1 



W-1 + 



w-1+1 

T E 



Q 



feo+3-1 



fcn = l 



Sm=l Qfco+m— 1 



(A.7) 



as well as the covariance 

Cov(pj,pi) = £[(pj - Mj)(p<- Mi)] 



1 / n 

^ k = l \%2m=l Qko+m-1 



Qk +i- 



k Sm=l Qk 



(A. 



+m-l 



In view of Eqs. (|A.2[) and (|A.3[ ). the quantities fij, Var(pj) and Cov(pj,p m ) are always finite irrespective of the existence 
of heavy tails in Qk which is for example the case of seismicity. Moreover, for the purpose of our calculations the 
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relation between the variance of pj, Var(pj-), and the covariance of pj and p m , Cov(pj,p m ), is important. Equations 
fO]) and CO)]) lead to 



(A.9) 



which when multiplied by (pj — fij) and averaged (cf. £ = w * l+1 Sfe =i ) results in 

Var(pj) = - ^ Gov(pj,p m ). 



(A.10) 



We now turn to the evaluation of the mean value /i of k\ obtained when the (natural time) window of length I 
slides through a time series of > 0, k = 1, 2, . . . W, 



W-l+l 



fj, = £(ki) 



W 



(A.11) 



k = l 



by studying its difference from the one that corresponds to the time series of the averages M = {nk} which is labelled 



i*i,m = E ( j ) Mi 



Ey^ 



(A.12) 



Hence, 



w-i+i i 



fi — K\M 



w 



S { E ^-bm(*o)-^ 



fco = l I m— 1 



-, 2 



E yP™( fc o) 



(A.13) 



In view of the definition of /Lt m , the first term in square brackets in the right hand side of Eq. (|A.13p vanishes, whereas 
the latter two terms reduce to the opposite of the variance of 



{x)m= E y^" 1 ' 



(A.14) 



leading to 



w-i+i ( i 



w - 



yyy E Ey \Pm(k ) - Mm] 



fe = l \m=l 

Expanding the term within the curly brackets and interchanging the summations, we get 

l 9 l-l I . 



Kl,M -M= E ~ Var (P™) + 2 E E IT Cov(pj,p m ). 
which, upon using Eq. (|A.10|) . leads to 



m— 1 j — 1 m— j+1 



M-Kl,M=X] E 



z / 



^Lcovfe !Pm ) = JEE 1Z -^Covfe, Pm ). 

3=1 m=3+l ^ 3=1 m=l 



The latter relation turns to 



E 

all pairs 



(j - m) 2 , . 

- Lov(Pj,p m ) 



(A.15) 



(A.16) 



(A.17) 



(A.18) 
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whci'C Eallpairs = E}=1 TL=j+1- 

Equation (|A.18[) shows that the mean value [i itself is a measure of the correlations between successive earthquake 
magnitudes. The practical use of this equation, however, in order to estimate the strength of these correlations between 
seismic events requires f32l l48M5l| the construction of a large number of shuffled copies of the original earthquake catalog 
and a comparison of (i with the relevant distribution obtained from the shuffled copies. Obviously, this task becomes 
cumbersome when the (natural time) window of length I is sliding through a long time series of Qk ■ 

When Qk are independent and identically distributed positive random variables, Eq. (|A.f 8|) results JlH, |49| in 



fj, = k u ( 1 - -J ) - k u (1 + I)Var(p), 



(A.19) 



where n u = 1/12 -corresponding to the K\ value for the uniform distribution- and Var(p) the variance of any pj. 



2. The standard deviation a of the m values 

Let us now investigate the standard deviation a of the n\ values obtained when the (natural time) window of length 
I slides through a time series of Qk- This is obtained from the variance 



u = Var (ki) = £ (k\ — m) 



1 



W-l+l 



W-1 + 



fcn=l 



(A.20) 



Numerically, the above quantity can be evaluated almost as easily as \i when Ki(fco) are available. 
In order to obtain an analytical expression, by inserting Eq. ([A~T8)) into (TA~20]) . we obtain 



w-i+i i 



W - 



7^1 E { E \^-\P™(ko)-H 
(j - mf 



ho — 1 \ m— 1 

2 



E -rP™( fc o) 



m— 1 
2 



E 

\m=l 



E 

all pairs 



Cov(pj,p m ) 



(A.21) 



Rearranging the terms 



we get 



X yP™( fc o) 



E 



\m—l 



E ~f [P™( fc o) - Mm] M E T bm(fco) + ^" 



= El bm(fco) - Mm] \ + 2{X)M \ E T bm(M ~ Vm] 



k rn—1 



(A.22) 



W-l+l 



W 



— r 

-i+i ^ 



fco=l 
\2 



1 ' 2 

m 



E ( JT ~ 2 (^T ) bm(fco) - Mm] ~ \ E T bm( fc 0) - Mm] 



E 1 S^ C 



all pairs 



(A.23) 
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Upon expanding the square over the square brackets in Eq. (|A.23[) we obtain six terms: 

2 



W 



1 W-l+l ( l , 2 \ "I ' 

yryr E |E (jr ~ 2 (x)mjJ b™0o) -Mm]j 



w-i+i r / 



W - 1 + 



E 

all pairs 



fco — 1 k. m— 1 
\2 



J ] X! ( ~ 2 (X>A4 — j k(fco) - Mm] 0]y bm(^o) ~ Mm] 



. m— 1 



(j - my 



Cov (pj,p m ) 



W -1 + 



W-l+l (l/o 

TTh 111 

-jY - 2(x)m— ) bm(fco) - 



fco=l ^ m— 1 



W 



~l + 1 E | E J ^(^0) - Mm] 



all pairs 



I 2 



-Cov(pj,p m ) 



W-l+l ( I 



W -1 + 



T E Ej b™( fc 0) - Mm] 



fc = l l.m=l 



E 

all pairs 



(j — J77,) 

^ Gov (Pj>Pm) 



(A.24a) 
(A.24b) 
(A.24c) 

(A.24d) 
(A.24e) 

(A.24f) 



The following comments are in order: First, the term in (|A.24cj) vanishes due to Eq. (|A.5[) . Second, the terms in 
()A.24bl) and (|A.24d[) clearly depend on moment correlations higher than the second, thus they should be neglected 
when restricting ourselves to second order correlations. Third, the second term in (|A.24eP can be evaluated using 
Eqs. (|A.15[) and (|A.17|) leading to a partial cancellation with the term in (|A.24fj) . Hence, restricting ourselves to second 
order correlations, we finally obtain 



W 



rf+T E |E (yr - 2 (^y) MM - m| 



E 



all pairs 



^ Cov(pj,p m ) 



(A.25) 



The first term in Eq. (|A.25[) can be evaluated by expanding the square over the curly brackets and using Eq. (|A.10p 
-in a way similar to Eqs. (|A.16p and (|A. 17|1 - so that we obtain 



a 2 = - 



E 

all pairs 



(7 "<*>«)' 



- ( y - (x)m 



Cov{pj,Pr, 



E 

all pairs 



(j — jjj \ 2 

^ Cov (PjiPm) 



(A.26) 



Equation (|A.26|) reveals that <r 2 -like the mean value in Eq. (|A.18p - is a measure of the correlations, but a 2 is almost 
proportional (see also below) to these correlations whereas in \i they appear as an additive term in Eq. (|A.18j) . 



3. The variability a//j, 

By combining Eqs. (|A.18j ) and (|A.26j) we find: 



E 



all pairs 



if -{x)m) 2 -{\-{x)m) 2 



COY(pj,Pr, 



E 



all pairs I 



(3 p Cov(pj,p m ) 



Kl,M + Eall pairs C ° V fe > P™) 



(A.27) 



This equation, which is just Eq.(4) of the main text, provides in general the interrelation between the variability (3 
and the event correlations. 

Additional insight on the physical meaning of a/ fx may be obtained when adopting the paradigm of the uniform 
distribution[32l. l36l |37| which corresponds to a simple system operating at stationarity, i.e., when Qk are independent 
and identically distributed positive random variables. In this case, we have[32j 



(A.28) 



9 



(x)m = 



mil 
~P = 2 + 21' 



(A.29) 



and due to Eq. ()A.10p 



thus wc obtain 



Cov(pj,p m ) = - 



Var(p) 



(A.30) 



a 2 = 



Var(p) 

(z=i) 



E 

all pairs _ 



m 1 1 

7 ~ 2 ~ 2/ 



2 2 22 



Var(p) 

(T=T) 



E 



(i - ™) 2 



all pairs 



r- 



2\ 



(A.31) 



For large I the summations over all pairs can be effectively, e.g. I > 10, approximated by integrations 



E 

all pairs 



mil 

7 ~ 2 ~ 27 



2 2 2? 



Z 2 ,i ,i 



^0 



■</• 



dxdip 



180' 



(A.32) 



E 

all pairs 



I 2 ~ 2 



o Jo 



12 



(A.33) 



Equation (|A.31|) simplifies to 



er 2 ps iVar(p)«^ 



- - ZVarfp) 
5 



and Ea. (|A.19|) becomes 

/i ps [1 - lVax(p)] . 
Thus, the variability simply results in 



/? = - = VN^pj 



| - War(p) 
1 - ZVar(p) 



(A.34) 



(A.35) 



(A.36) 



When Qk exhibit heavy tails as in the case for seismic- 
ity, the quantity lVai(p) measures the intensity of such 
tails and so does (3. In other words, for randomly shuf- 
fled earthquake data or earthquakes occurring with tem- 
porally uncorrelated magnitudes, the variability /3 is a 
measure of the 6-valuc of the Gutcnbcrg-Richtcr law. 
Since real seismic data may also exhibit temporal correla- 
tions between earthquake magnitudes (27l. l32l I48l - l52j . the 
general expression of the variability obtained above (i.e., 
Eq. (|A~27)) ) from Eqs. (jA~l~8|) and (|A~26f capt ures both the 
effects of correlations and heavy-tails. 



When do not exhibit heavy tails, which is not of 
course the case of seismicity, the quantity lV&r(p) is sim- 
ply related [12, HH to the mean fj. and the standard de- 
viation ctq of Qk'. 



(A.37) 



«Var(p) = H. 

Assuming that oq/hq is of the order of unity, ZVar(p) 
becomes small compared to unity when I > 10, and 
Eq. (|A.36[) becomes 



B = - = —— (— 

M \/5Mo \y/l 



(A.38) 



i.e., the variability of «i is directly proportional to the 
variability of the data Qfc- Note that the same holds for 
the standard deviation of the natural time entropy [321 , |53| 
S as well as for change AS of the entropy in natural 
time under time reversal [32l (cf. for the analysis in 
natural time under time reversal, see also Refs.[55[ and 
[56|). Thus, one could alternatively view (3 as an entropic 
measure. 
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